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ABSTRACT 

A hyperaccretion disk formed around a stellar mass black hole is a plausi- 
ble model for the central engine that powers gamma-ray bursts (GRBs). If the 
central black hole rotates and a poloidal magnetic field threads its horizon, a 
powerful relativistic jet may be driven by a process resembling the Blandford- 
Znajek mechanism. We estimate the luminosity of such a jet assuming that 
the poloidal magnetic field strength is comparable to the inner accretion disk 
pressure. We show that the jet efficiency attains its maximal value when the 
accretion flow is cooled via optically-thin neutrino emission. The jet luminosity 
is much larger than the energy deposition through neutrino-antineutrino anni- 
hilation [yv — > e + e~) provided that the black hole is spinning rapidly enough. 
When the accretion rate onto a rapidly spinning black hole is large enough 
(> 0.003 — O.OlM s _1 ), the predicted jet luminosity is sufficient to drive a GRB. 

Subject headings: accretion, accretion disks - black hole physics - gamma rays: 
bursts - neutrinos 

1. Introduction 

Gamma-ray bursts (GRBs) are the most luminous objects in the Universe, releasing 
> 10 51 erg in a few seconds. The spectral features and lightcurves of their prompt and 
afterglow emission imply that they are produced in ultrarelativ istic jets. The central engines 



most likely to launch these jets are hy per accreting black holes ( INarayan. Paczyhski fc Piran 



19921 ; INarayan. Piran fc Kumarll200ll ). In this model, a relativistic jet emerges from a system 



with a massive accretion disk (0.1 — 1M ) surrounding a stellar-mass black hole, a system 
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expected to form after a cataclysmic event such as th e gravitation al collapse of a massive 
star or a merger of a neutron star-neutron star binary (|Piranlll999l ). 



The observed luminosity requires a very large mass accretion rate, ~ 0.01 — IOMqS -1 . 
In such a case the accretion flow is extremely optically thick with respect to photons and it 
cannot cool efficiently via radiation. However, since the density and temperature are very 
large (p > 10 7 g cm" 3 , T > 10 10 K) when the mass accretion rate is sufficiently high, neu- 
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The accreting system must launch a relat i vistic jet to produce a GRB. Neutrino pair 
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Nagataki 2009) are the two 


processes most discussed in the literature. The goal of this work 



is to estimate the luminosity of the jets emerging from a hyperaccretion flow driven by an 
MHD mechanism, as well as the luminosity of thermal neutrinos emitted from the disk. The 
first is estimated using the maximal magnetic field sustainable on the event horizon, which 
is expected to be limited by the pressur e in the disk near the ISCO, the innermost stable cir- 



cular orbit (IKrolik fc Piranl l201ll |2012| ) . The second we calculate directly from the density, 



pressure, and lepton number of the material in the disk. 

In Section 2 we predict the jet power in the MHD model by solving for the structure of a 
surrounding NDAF and compare this luminosity to the energy deposition rate via neutrino 
pair annihilation above the accretion flow. We discuss the parameter dependence of our 
model and its application to GRB jets in Section 3. We summarize our results in Section 4. 



Model 



We consid er a rotating black hole wh ose horizon is threaded by a large-scale poloidal 
magnetic field ( IBlandford fc Znajeklll977l ). The Poynting luminosity expected from such a 



system is ~ c(B 2 /8ir)R^, where B is the poloidal magnetic field strengt h on the horizon 
and R g = GM-q-r/c 2 is the gravitational radius of the central black hole ( iMcKinneyl 12006 
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Hawley &: Krolikl 120061 ). Following iKrolik fc Piranl ( 1201 ll ) we estimate the jet luminosity as 



f(a/M BH )c(B 2 /8n)R. 



(1) 



where /(o/Mbh) is an incre asing function of \a/M^\ whose exact form depends on the spe- 



cific field configuration (e.g. I Blandf ord fc Znaiek 



20081 ; iTchekhovskoy. Narayan fc McKinneyl l201ll ) . The re have been some attempts to calcu 



19771 ; iTchekhovskoy. McKinney fc Narayan 



late / "(q/Mrh) analytically for specified field geometries (IBlandford &; Znajek!ll977l : ITchekhovskoy. Narayan 
20101 ) . but these configurations were not required to be consistent with the dynamics of the 
surrounding accretion flow. They were also relatively simple poloidal topologies, but as 
Beckwith. Hawley &; Kroliki (120081 ) showed, high jet power requires that the vertical compo- 
nent of the magnetic field maintain the same sign for long times; in other words, low-order 
poloidal topologies are the most favorable to strong jets. Consequently, it is possible for 
/(o/Mbh) to vary over an extremely wide range, from very small numbers to ~ 0(1). When 
the topology is optimal and a/M > 0.9, the /(o/Mbh) arising from a dynamically self- 
consistent field structur e can be > 0.05, and even larger for spin parameters closer to unity 
rtrlawlev fc Kroliki l2006h . 



Beckwith. Hawley fc Kroliki ( 120091 ) argued that the magnetic pressure near the horizon 



may be limited by the inner disk pressure. The magnetic field energy may therefore be 
estimated (or at least bounded) by the disk pressure pdisk near the ISCO: 



8tt 



Pdisk(-Ri) 



(2) 



where i? in is the radius of the pressure reference point and is the ratio of the midplane 
pressure at R{ n to the magnetic pressure in the black hole's stretched horizon. 

To evaluate the midplane pressure of the hyperaccretion flow, we solve for the disk 
structure in the innermost region (R ~ 0(R g )), where the jet is expected to be launched. 
As we show below, the behavior of an hyperaccretion disk is mainly determined by its mass 
accretion rate M. To highlight what we believe to be the principal physical mechanisms 
(and also eliminate mathematical complications), we adopt a highly simplified model for 
this structure in the expectation that the trends we uncover while exploring a factor of 10 5 
range in M will be robust with respect to refinements in the model. In particular, we assume 
that both the disk dynamics and the gravitational potential are Newtonian despite the fact 
that this near a black hole general relativistic effects are important. We also assume that the 
system is stationary and axisymmetric. We adopt the Shakura-Sunyaev formalism for the 
disk structure, which formally applies only to thin disks, assumes that the inflow speed is 
very small compared to the orbital speed, and parameterizes the r-0 component of the stress 
( IShakura fc Sunyaevlll973l ). When the disk becomes geometrically thick, in the advection- 
dominated regimes, we use in our analytic estimates (but not in the numerical solution) a 
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dimensional analysis. Remarkably the results are similar up to a numerical coefficient to 
those obtained using the Shakura & Sunyaev model. We further set the correction factor X 
appearing in Equation H] to always be unity, when it should in fact be somewhat less than 
that. These are rough approximations, but we nonetheless expect our results to give the 
right qualitative picture in this scenario. 

Here we present the basic disk equations for the density p, the temperature T, and 
the scale height H. These are the expressions for mass conservation, angular momentum 
conservation, energy conservation, and hydrostatic balance applied at i? in : 

M = -2nR m Zv R , (3) 
2aHp Aisk = Mn } Rm) X(R/R g ), (4) 

Q + = Q~, (5) 



Pdisk 
P 



n{R hi ) 2 H\ (6) 

(7) 



where S, vr, Q(R), and a denote the surface density (= 2pH), radial velocity, angular 
velocity (= GM Bii /R 3 ), and ratio of integrated stress to integrated pressure, respectively. 
In Newtonian dynamics, the factor X(R/R g ) accounts for the reduction in stress due to the 
net angular momentum flux through the disk; although the net angular momentum flux per 
unit mass accreted has traditionally been iden tified with the specific angular momen tum of 
orbits at the ISCO INovikov fc Thome! Jl973|) , it a lso depends upon MHD effects (IKrolikI 



1999; Gammie 1999 



Krolik. Hawlev fc Hirose!l2005h . 



The pressure pdisk is the sum of the contributions from radiation, baryonic gas, degen- 
erate electrons, and neutrinos (if they are trapped): 



4/3 

4/3 



11 r^A P^bT 2ixhc ( 3 , 
Pdisk ^ -aT 4 + ^-^ + — (Y e p) 

12 m p 3 \QTimpJ 

where a ~ 7.56 x 10 _15 erg cm -3 deg -4 is the radiation constant, kg is the Boltzman constant, 
and m p is the mass of a proton. The first term on the right hand side represents the 
contribution from both photons and relativistic electron-positron pairs; the coefficient is 
11/12 when the temperature is higher than ~ 10 10 K. In the third term, which corresponds 
to degeneracy pressure, Y e is the electron fraction, which should be evaluated from the 



-5 - 



/3-equilibrium condition0 



l-Y e 



exp 



Q-fJ- t 
' k B T 



(9) 



where Q ~ 1.29MeV is the difference between the neutron and proton masses, and fi e is 
the chemical potential of the electrons (we assume that the chemical potential of neutrinos 
can be ignored an d that heavy nuclei are totall y dissociated in the innermost region of the 
accretion flow; see lKawanaka fc Mineshigdl20071 ). The last term represents trapped neutrino 
pressure; u v is the neutrino energy density, defined below. 

Q + is the heating rate per unit area: 



cr 



3GM RH M 



AttR 3 

9.8 x 10 36 erg cm" 2 s _1 (M/Kr 3 M s" 
x(M BH /3M )- 2 (i?/6iy- 3 . 



(10) 



Q~ is the cooling rate per unit area, including both neutrino and advective cooling: Q~ = 
Q~ + Q~ dv - The advective cooling rate Q~ dv , which equals the inward flux of the gas energy 
along the flow, can be described as: 



TY.VR-. 

-Kin 



'IV. 



Here s denotes the entropy per unit mass. There are se veral significant neutrino emis sion 
processes expected in such a hyperaccretion flow (see also iKawanaka fc Mineshigd 120071 and 
references therein): electron/positron capture on nuclei (p + e~ — > n + v e \ n + e + — > p + z/ e , 
also known as the Urea process), electron-positron pair annihilation (e~ + e + — > vi + 
where i represents electron-, mu-, and tau-type neutrinos), nucleon-nucleon bremsstrahlung 
(n + n — > n + n + Ui + Ui), and plasmon decay (7 — > i/ e + u e ). We should also take into 
account the neutrino optical depth in the accretion flow, which can be larger than unity 
when the mass accretion rate is sufficiently large. Here we adopt a prescription based on 
the two-stream approximation (jPi Matteo. Perna fc Narayanll2002l ; iKohri. Narayan &: Piran 



1 Note that this approximation is valid only when the flow is optically thick with respect of v e , and 
when neutrino emission is not efficient the electron fraction should be nearly equal to 0.5. In the numerical 
calculations we make use of a bridging formula in evaluating Y e to connect these two extreme regimes 
smoothly. 
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20051 ): 

Q- = 2 V {V8)aTA (12) 

Y(3/4)(W2 + l/^+l/3rJ' 

= E W + i/^) 

Y W2 + l/V3+l/3 Vi 

where = r a ^. + r S)1/j and the neutrino flavor is labeled by subscript i = (e, fi, r). T a>v . and 
T s ,ui represent absorptive and scattering opt ical depths respectively. The d etailed description 



of these neutrino opacities can be found in iKawanaka fc Mineshigd (120071 ). 



When the optical depth is so large that the neutrino diffusion timescale in the accretion 
flow is longer than the accretion timescale, neutrinos are completely trapped in the accretion 
flow. In that state, advective cooling dominates: Q~ ~ Q&dv- The accretion flow would also 
be advection-dominated when the mass accretion rate is small enough that the density and 
temperature of the disk are so low that neutrinos are not efficiently emitted, but is still much 
larger than the Eddington accretion rate. 

We can divide the hyperaccretion flow into several regimes according to the dominant 
physical processes taking place. In particular, the cooling processes (advection, optically- 
thin neutrino emission and optically-thick neutrino emission) and the pressure components 
(radiation, baryonic gas, and degenerate electrons) vary from one regime to another. In 
the rest of this section, we solve the structure of a hyperaccretion flow in its innermost 
region as a function of the mass accretion rate and examine the dominant physical processes 
determining the behavior of the accretion flow in the different regimes. 



2.1. Numerical Results 

We numerically solve the fundamental equations of an hyperaccretion flow, Eqs. j3J)-([S]), 
in the innermost region to find the density, temperature, and scale height as functions of the 
mass accretion rate M. In the course of doing so, we also identify the dominant physical 
mechanisms operating over a wide range of M values. As the relevant ranges of parameters 
for the central engine of a GRB, we consider 0.01 < a < 0.3, 1 < Mbh/Mq < 30, and 
2 < Ri n / R g < 6. Our fiducial parameters are R in = 6R g , M B h = 3M and a = 0.1. 

Figure [TJ depicts the different cooling rates (neutrino emission and advection), optical 
depths for different flavors of neutrinos, and pressure components (radiation, baryonic gas, 
and degenerate electrons) as functions of mass accretion rate. We can see that neutrino 
cooling becomes dominant over advective cooling above M ~ O.OO3M s _1 . 
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mass accretion rate [M s" ] 

Fig. 1. — Physical properties in the innermost region of a hyperaccretion flow as functions of 
the mass accretion rate. Top: cooling rates by advection (solid line) and neutrino emission 
(dashed line). Middle: the optical depths for electron neutrinos (thick solid line) and fx/r 
neutrinos (dashed line). The boundary between optically thin and thick is also shown (thin 
solid line). Bottom: the pressure from baryonic gas (solid line), relativistic particles (photons, 
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Fig. 2. — Jet luminosities and efficiencies as functions of mass accretion rate calculated 
numerically (solid line) and analytically (dashed line). Here (3h and /(o/Mbh) are set to 
unity. 
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Figure |2] depicts the jet luminosity L )et and efficiency rj- ]Ct = Lj et /Mc 2 , estimated from 
Eq. ([1]). Inspection of Figures [1] and [2] reveals a sharp break in the nature of the disk and jet 
at M ~ 0.003M Q s _1 , the transition from advection to neutrino cooling. In order to interpret 
this break, we have constructed an approximate analytical model for hyperaccretion disk 
structure. This model can qualitatively reproduce not only the physical properties of the 
flow, but also the mass accretion rates at which one regime changes to another. 



2.2. The Analytic Model 

2.2.1. p ~ p ra( j and advection cooling (M < O.O18M s _1 for our fiducial parameters) 

When the mass accretion rate is not large enough for the disk to cool efficiently via 
neutrino emission, the accretion flow is advection-dominated, and the pressure is dominated 
by photons and relativistic electron/positron pairs. In this regime, the accretion flow is 
geometrically thick and the thin-disk approximation is not valid. Instead, we can estimate 
the disk pressure using dimensional analysis^]. Because the aspect ratio of the accretion flow 
H/R ~ 1, the density at R in is ~ M / (27iRf n Vfi). As the gas does not cool efficiently, the ratio 
of pressure to density (i.e. the square of the sound speed) is close to the virial temperature 
ksTvir ~ GM-Biimp/ Ri n . Using the relation between vr and the circular orbital speed t> or b, 
Vr/v ox \> ~ a, the pressure can be evaluated as: 

Pdisk(-Rin) ~ PC 2 S 

M / GMbh V 72 
2nRla \ R iQ J 
~ 5.4 x 10 27 dyn cm" 2 o^j(M /lO" 3 M s _1 ) 

x(M BH /3M )- 2 ( J R in /6 J R 9 )- 5/2 . (14) 

Using this expression for p^isk we estimate the jet luminosity as 

L jet ~ 3.2 x 10 49 erg s" 1 ^^(M/IO^MqS- 1 ) 

x(i? in /6/y- 5 / 2 /(a/M BH ). (15) 

The corresponding energy efficiency of the jet: 

~ 0.018/3 h - 1 a ;}( J R in /6i? ;? )- 5/2 /(a/MBH) (16) 



2 Remarkably, the results are qualitatively similar to those obtained using the thin disk approximation. 
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is independent of the mass accretion rate. 

The accretion flow is advection-dominated as long as the neutrino emissivity per unit 
volume is smaller than half the advective cooling rate. In this regime, neutrino emission is 
dominated by electron/positron capture, and the neutrino emissivity q~ is given by 

q~ ~ 9.0 x 10 33 erg cnT 2 s p 10 T^, (17) 

where p w = p/10 10 g cm" 3 and T n = T/10 n K (Popham, Woosley & Fryer 1999). The 
density and temperature of the gas are: 



P 



M 



2nR( n av orh 
~ 3.6 x 10 7 g cm" 3 a^^M/lO^MoS" 1 ) 

x(M BH /3M )- 2 ( J R in /6 J R,)- 3 / 2 , (18) 

i2 Pdisk y /4 

11a J 

~ 3.0 x 10 10 K ^^(M/IO^MqS- 1 ) 1 / 4 

x(M BH /3M )- 1 / 2 ( J R in /6i? s )- 5/8 . (19) 

Substituting these values into the expression for neutrino emissivity, we obtain: 

q~ ~2.3 xl0 28 erg cm" 2 s' 1 a^ /2 (M/lO- 3 M s _1 ) 5/2 

x (M BH /3M G y 5 (R in /QR g )- 21 /\ (20) 

On the other hand, the advective cooling rate per unit volume is 

^ v - Tv RP ^ 

~ 3.6 x 10 30 erg cm' 3 s" 1 (M/IO^MqS^ 1 ) 

x(M BU /3M Q )' 3 (R in /QR g )-\ (21) 

Here we use s = 'jpdisk/ pT(j — 1) where 7 = 4/3 is the adiabatic index of relativistic particles. 

A transition takes place when q~ > (l/2)q~ dv , i.e. 

M/Mqs" 1 ~ O.O18^ / 1 3 (M BH /3M ) 4 / 3 ( J R ill /6i?,) 5/6 . (22) 

Above this accretion rate the accretion flow is no longer advection-dominated, as neutrino 
cooling determines the structure of the flow. In addition, the dominant pressure source is 
no longer radiation, but baryonic gas pressure. In the following, we call this transition mass 
accretion rate M ign ("ign" stands for ignition of neutrino emission). 



- 11 - 



2.2.2. p ~ Pgas and optically thin v e cooling (O.O18M s 1 < M < O.O45M s 1 /or our 

fiducial parameters) 

In this regime the disk is mainly cooled by neutrino e mission via electron /positron 



capture, while the pressure is dominated by baryonic gas (see iKawanaka &: Mineshigd 12007 



and references therein). In this case, the pressure can be expressed as 

M* - (23) 
m p 

The neutrino cooling rate per unit surface area by electron/positron capture is given by 
Eq. ( TT7j) . Substituting Q~ = Q~ in Eq. ([5]), we solve the basic equations and derive the disk 



properties in this regime: 



p = 6.3 x 10 10 g cm" 3 a'f^M/O.OlMoS" 1 ) 

x (Mbh/SMo)- 17 / 1 ^^^^)" 51 / 20 , (24) 

T = 3.6 x 10 10 K al^(M BU /3M Q )- 1/5 

x(R in /QR g )- 3 / 10 , (25) 

H = 3.8 x 10 5 cm al[l\M Bn /3M Q f /w 

x(R in /QR g ) 27 / 20 . (26) 

The disk pressure is given by 

p disk ~ 1.9 x 10 29 dyn cm -2 a ai 11/lo (M/O.OlM s _1 ) 

x(M BH /3M )- 19 / lo (i? in /6 J R ;? )- 57 / 20 . (27) 

Finally we estimate the jet luminosity: 

L jet ~ 1.1 x 10 51 erg s" 1 ^^^(M /OMM^ 1 ) 

x(M B H/3M ) 1 / lo ( J R in /6i? 9 )- 57 / 2O /(a/M B H). (28) 

The jet luminosity in this regime depends linearly on the mass accretion rate. The efficiency 
7/ jct = L je t/Mc 2 is: 

r^VlO/^/o^U/lO/p. , R „ X-57/20 



77jet ~ 0.061^^-1 /lu (M BH /3M ) i / iU ( J R in /6i? 



9 J 



x/(a/M BH ), (29) 

and is independent of the mass accretion rate. The higher jet luminosity and efficiency in 
this regime is the result of more efficient neutrino cooling, which leads to a geometrically 
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thinner accretion flow. Thinner disks accrete more slowly. In the end, the midplane pressure 
Pdisk increases, permitting it to confine a stronger magnetic field on the black hole's event 
horizon. This stronger magnetic field supports a more powerful jet. 



The system is very optically thin to neutrinos at this stage: 

r a ,, e ~ 4.5 x 10- 7 T? lPw H, 



2.7 x IQ- 7 T? lPl0 H 



(30) 
(31) 



where r a ,, y and the optical depth due to nucleon absorption and scattering, respec- 

tively (see lDi Matteo. Perna &: Narayanll2002l ) . The total optical depth for electron neutrinos 
T Ve = r a ^ + T SVe exceeds unity only when 



M/Mos" 1 > O.O45ao / i 5 (M BH /3M ) 6 / 5 
x(R m /QR g ) 9 /\ 



(32) 



and then our expression for Q~~ would not be valid. This mass accretion rate marks the 
transition to the third regime, which is optically thick to v e . In the following we call this 
transition M t hi C k- 



2.2.3. p ~ Pgas and optically thick v e cooling (t).O45M s 1 < M < 2.2M Q s 1 for our 

fiducial parameters) 

In this regime the innermost region of the disk is optically thick to electron neutrinos, 
and the opacity is dominated by absorption onto nucleons. While the disk is not so optically 
thick that the neutrinos are totally trapped in the flow, the cooling rate is: 

QZ ^ ^-^1, (33) 

where t v is the neutrino optical depth. 

First we consider the case in which the disk is optically thick only to v e and thin for 
and v T . With cooling only by electron neutrino emission, the cooling rate is 

QZ ~ 7.4 x 10 34 erg cm" 2 s' 1 T^p^H' 1 , (34) 

which should be equated with the heating rate in the energy balance equation. Solving the 
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fundamental equations of the disk we derive its density, temperature and scale height: 

p « 4.1 x 10 n g cm" 3 a~i /2 (MBH/3M )- 1/2 

x(R in /6R g )- 3 / 4 , (35) 
T « 4.8 x 10 10 K a^ 1 l /3 (M/0.1M Q s- 1 ) 2/3 (M B}1 /3M & )- 1 

x(i? in /6i? 9 )" 3/2 , (36) 
H « 4.4 x 10 5 cm (^/OTM s-y /3 (M BH /3M ) 1/2 

x(R in /6R g ) 3 / 4 . (37) 

The disk pressure in the midplane is 

Pdisk » 1.6 x lO 3O a o ^ /6 (M/O.lM s- 1 ) 2/3 (M BH /3M )- 3 / 2 

x( J R in /6i? 9 )" 9 / 4 , (38) 

and the jet luminosity is 

L jet ~ 9.8 x 10 51 erg s" 1 /3^ 1 a o ^ /6 (M/O.lM s- 1 ) 2/3 

x(M BH /3M ) 1 / 2 ( J R in /6 J R 9 )- 9 / 8 /(a/MBH). (39) 

The dependence on the mass accretion rate is weaker than in the last regime. The efficiency 
is: 

%t ~ O.O55^ 1 a o ^ /6 (M/O.lM s- 1 )- 1/3 (M BH /3M ) 1 / 2 

x(R in /QR g )- 9 ^f(a/M BU ). (40) 

In order for the disk to cool by neutrinos efficiently, the neutrino diffusion time tdifl> e = 
Hr Ue /c must be shorter than the accretion time t acc = 1/ '(ctQ)(R/ ' H) 2 . This condition can 
be expressed as 

M/MqS" 1 < 2.2^/ 1 16 (M BH /3M ) 21 / 16 

x(i? in /6i?,) 51 / 32 , (41) 

and when the mass accretion rate exceeds this value, electron neutrinos are no longer able 
to escape the disk. We call it M tra _ VyUc . For the other kinds of neutrinos (z/ M and u T ), the 
diffusion time is still shorter than the accretion time, so they can cool the disk efficiently. 
This is the fourth regime, which we discuss in the next subsection. 



2.2.4- P ~ Pgas, and optically thick v x cooling (2.2M Q s 1 < M < 4.1M s 1 for our fiducial 

parameters) 

Even when the mass accretion rate is so large that z/ e s are completely trapped in the 
disk, z/^s and u T s emitted via the electron-positron pair annihilation process may escape the 



-14- 



disk so long as their diffusion time is shorter than the accretion time. For those neutrinos, 
the opacity is dominated by scattering on nucleons (Eq. l3Tj) . Then the cooling rate is 

« , 2x (2 4{7/2) ° Ti 



~ 2.0 x f0 35 erg cm" 2 s" 1 T^p^H, (42) 

where the factor 2 in the first line represents the number of flavors of neutrinos that we 
are interested in {v^/v^ and u T /u T ). Substituting this Q~ into Eq. (jSJ), we derive the disk 
properties: 

p ~ 9.3 x 10 u g cm" 3 o;^i /2 (Mbh/3M )- 1/2 

x(R in /6R g )~ 3 /\ (43) 
T ~ 1.3 x 10 U K a a i /3 (M/lM s~ 1 ) 2/3 (M BH /3M )- 1 

x(R in /QR g )~^ 2 , (44) 
H ~ 7.3 x 10 5 cm a a i /6 (M/lM s" 1 ) 1/3 (M BH /3M ) 1/2 

x(i? in /6i? 9 ) 3 / 4 . (45) 

The disk pressure is then 

p disk ~ 9.9 x 10 30 dyn cm" 2 a^ /6 (M/lAf Q s -1 ) 2/3 

x(M BH /3M )- 3 / 2 ( J R in /6i? ;? )- 9 / 4 , (46) 

and the jet luminosity is 

L jet ~ 6.0 x 10 52 erg s" 1 ^^^(M/lAfoS" 1 ) 2 / 3 

x(M BH /3M ) 1 / 2 ( J R in /6i? ;? )- 9/4 /(a/MBH). (47) 

The luminosity now depends on the mass accretion rate in the same way as in the last regime. 
The efficiency is 

Viet ~ O.O33/3 A 1 « o ^ /6 (M/lM s- 1 )- 1/3 (M BH /3M ) 1 / 2 

x(R in /6R g )~ 9 ^f(a/M BR ). (48) 

The timescale for p/r neutrinos to escape the disk diffusively is 

Ht 8 „ 



t. 



diffj 



C 

5.2 x 10" 4 s a^ /2 (M/lM s- 1 ) 2 (M BH /3M )- 3/2 
x(iW6iy- 9 / 4 , (49) 
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and the accretion timescale is 

1 



^acc — 



R^ 2 



<e) 2 



~ 3.9 x l(T 2 s a^ /3 (M/lM B- 1 )- 2/3 (M BH /3M 

x(R in /6R g ) 2 . (50) 

Therefore, fi/r neutrinos can escape an accretion flow without being trapped as long as 

M/MqS- 1 < 4.1^/ 1 16 (M BH /3M ) 21 / 16 

x(R in /QR g ) 51 / 32 . (51) 

At this mass accretion rate, there is a transition to the fifth and last regime, in which the 
accretion flow is once more advection-dominated. 



2.2.5. p ~ p ra d +p e ± + ^Pua Q ~ Qa,dv f4.1M s 1 < M for our fiducial parameters) 

When the mass accretion rate is so large that the all kinds of neutrinos are trapped, the 
accretion flow cannot cool and the flow is advection-dominated. In this case the pressure 
would be dominated by several kinds of relativistic particles, photons, relativistic electron- 
positron pairs, and all three flavors of neutrinos. We can estimate the disk pressure in a way 
similar to Sec. 2.2.1 to find 

p disk ~ 2.7 x 10 31 dyn cm" 2 a^(M/lOAf s -1 ) 

x(M BH /3M )- 2 (iW6iy- 5 / 2 . (52) 

The jet luminosity is then 

L jet ~ 3.2 x 10 53 erg s" 1 ^ l a^(M /lOA^s" 1 ) 

x(R m /QR g )- 5 / 2 f(a/M B]1 ). (53) 

As a result, the efficiency is the same as in Equation ( 1161) . 

r? jet ~ 0.018/3- 1 a :}( J R in /6i? g )- 5 / 2 /(a/M BH ). (54) 



2.3. The Order of the Transition Points 



In the discussion above, we have picked out the four mass accretion rates that divide 
the five states of an hyperaccretion flow. However, the sequence we have presented does not 



0.01 \r 

0.001 



ADAF to thin NDAF 
thin to thick NDAF 
v e trapped 
v x trapped 



_i i i i i i i_ 



1 



10 

black hole mass [M Q ] 



Fig. 3. — The four transition mass accretion rates as functions of the black hole mass. Here 
we fix the other parameters as a = 0.3 and R[ n = 6R g . Note that, for canonical values of 
these parameters (i.e. a = 0.1 and Ri n = 6R g ), there is no change of the order of transition 
points in this black hole mass range. 

always apply. These mass accretion rates are functions of a, Mbh, and R{ n , and in some 
regions of parameter space, the transition from advection-dominated flow to neutrino-cooled 
flow (Mjga) takes place at a higher accretion rate than the transition to optically thick flow 
(with respect to u e ). This means that the second regime of optically thick neutrino cooling 
does not exist for these parameters. For example, we show the transition mass accretion 
rates as functions of the black hole mass (while the other parameters are fixed at a = 0.3 
and -Ri n = QR g ) in Figure |3j One can see that there is a critical black hole mass at which 
the first transition and the second transition change places, and also that for this choice of 
parameters the mass accretion rate range of the optically-thin NDAF regime is quite narrow, 
if it exists at all. In Figure H^ve show Mthick/-^ign — 1 f° r various parameter sets as functions 
of the black hole mass. For (R in ,a) = (GR g , 0.1) and (4R 9 ,0.1), the curve is well above 
^thick/A^ign — 1, and the optically-thin NDAF regime exists for each of those parameter sets. 
The general condition for such a change in the order of the transition points is: 

a 13/15 (M BH /M Q ) 2/15 (R/R g )- 29/30 > 0.07. (55) 

When this condition is satisfied, advection-dominated flow transforms directly to optically 
thick neutrino cooling at ( 132]) . Figure (5Uisplays the limited range of parameters satisfying 
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Fig. 4. — M t h ic k/M ign — 1 as a function of the black hole mass. Black and red lines represent 
the case with a = 0.1 and 0.3, respectively, and solid, dashed and dotted lines represent the 
case with R[ n = 6R g , 4R g and 2R g , respectively. The horizontal thin black line represents 

Mthick/Mgn = 1. 



this condition. In other words, except in the case of extreme values of the parameters, the 
five states for hyper accretion flows are ordered as we have described them. 



3. Results and Discussion 

3.1. Parameter Dependence of the Jet Luminosity 

From the analytic model presented above, we can estimate the jet luminosity and ef- 
ficiency expected from a hyperaccreton flow as functions of mass accretion rate for various 
choices of parameters (a, M B h, Rin)- Figure |6] shows how the jet luminosity based on the BZ 
scenario (thick lines) depends on the mass accretion rate. The different lines show the results 
for different parameters. We can see that with all parameter sets, the jet luminosity has a 
significant discontinuity at M; gn . This jump, as explained already, arises from the transition 
between an advection-dominated accretion flow, in which neutrino cooling is not efficient, 
and a neutrino-dominated accretion flow. In the top panel, we can see that for larger black 
hole mass the mass accretion rate at which the transition from the advection-dominated flow 
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Fig. 5. — The optically thin neutrino cooling regime exists above (and to the right) of the 
marked plane in the parameter phase space of (a, Mbh, Rm)- 

to the neutrino-cooled flow (M; gn ) is larger. The value of M- lgn rises from ~5x 10 _3 M Q s _1 
at M B h = 1M© to ~ O.3M s _1 at M B h = 3OM . This is because the temperature of the 
accretion flow is lower for larger black hole mass when the radius (in units of the gravita- 
tional radius) and the mass accretion rate are fixed. Therefore, a larger mass accretion rate 
is needed for an accretion flow around a larger black hole to cool efficiently via neutrino 
emission. 

A similar explanation also applies to the bottom panel, which shows the dependence 
of the jet luminosity on R in . Here, for smaller i? in the temperature of the accretion flow is 
higher, which results in a smaller Mi gn and larger luminosity. Over the range of plausible 
values of R[ Q , the variation of M ign is much smaller than for the dependence on Mbh'- it 
changes by only a factor ~ 4. 

The middle panel shows that the jet luminosity depends quite strongly on the stress/pressure 
ratio in the sense that for larger a, the luminosity is smaller; equivalently, M ign increases with 
increasing a. Over a span of a factor of 30 in a, from a = 0.01 to a = 0.3, M ign increases 
by a factor of 300. To interpret this behavior, we note that the inflow rate is proportional 
to a (i.e. the accretion timescale is proportional to a -1 ). Therefore, when a is larger in 
the radiation-dominated advective regime, the midplane pressure is smaller (eqn. [T4"|) . The 
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0.0001 0.001 0.01 0.1 1 10 

mass accretion rate [M s" 1 ] 



Fig. 6. — Jet luminosity expected from the BZ mechanism (thick lines) and neutrino pair 
annihilation (thin lines). Their dependence on the central black hole mass {top), on a 
(middle) and on the inner radius of an accretion flow (bottom) are shown. The model 
parameters are a = 0.1, M B h = 3M Q and R m = QR g unless otherwise stated. Here /3h and 
f(a/M BH ) are set to unity. 
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temperature is also lower (eqn. 120]) . These two trends together diminish the neutrino cooling 
rate, forcing the transition accretion rate higher. They are reinforced by the scalings on the 
opposite side of the transition, in which p oc a _13//10 , while the temperature depends much 
more weakly on a (only oc a 1 / 5 ). Thus, the jet luminosity is strongly reduced by a larger 
stress/pressure ratio. 

Figure [7] similarly shows how the jet efficiency depends on our three free parameters. 
As we have previously emphasized, the principal dependence is on accretion regime: the jet 
efficiency is greatest when the flow is in the optically thin neutrino-cooling state. Although 
this state occurs over different ranges of accretion rate for different black hole masses, the 
efficiency achieved there varies rather little with Mbh- On the other hand, even after allowing 
for the different accretion rate ranges in which the flow is neutrino-cooled, it is quite sensitive 
to both a and R in / R g . In that favored regime, ?7j et is roughly oc a~ l . Still more dramatically, 
a factor of 3 decrease in R in / R g leads to a factor of 30 increase in rj )et at the respective peaks 
of jet efficiency in the optically-thin neutrino-cooled state — and this comparison is taken 
at fixed /(o/Mbh), whereas realistically one would expect smaller R m /R g to be associated 
with more rapid spin and therefore a significant increase in /(o/Mbh)- Note that in a 
small region of the lower panel r]j et > 1. While possible, in principle, in a BZ process, our 
Newtonian approximation is not suitable to explore this regime. 



3.2. BZ versus Neutrino Annihilation Jets 

Neutrino-antineutrino annihilation i nto electrons and p ositrons (hereafter "neutrino pair 



annihilation") above the accretion flow (jEichler et al.lll989l ) is another mechanism proposed 



for the energy source of relativistic jets that power GRBs. It is worth comparing the lumi- 
nosity and efficiency of such neutrino pair annihilation driven jets with the BZ driven jets 
expected from our NDAF model. 

In estimating the energy deposition rate due to neutrino pair annihilation above a 
hyperaccretion flow, we adopt a simplified formulation: the accretion flow is assumed to be 
geometrically thin (i.e. neutrinos are emitted from the equatorial plane), and we neglect the 
bending of neutrino paths and the gravitational redshift. The energy deposition rate per 
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unit volume above the accretion flow can then be calculated from 

AL k AL-' 



fe "fc fc , a fc' \ e M\ e M 



(56) 



where ^ = ( 7 (C 1 +C , 2 ) i , i p i /(127r 2 c 5 m 2 ) « 1.7x 10- 44 cm ergs" 2 s" 1 and A 2 = <7 Q (C 3 W( 4 ^ 2 c) 
1.6 x 10~ 56 cm s _1 are the neutrino cross section constants for electron neutrinos (Ruffert 
et al. 1997). Here AL^. and (e) Ui are the neutrino luminosity coming from the fc-th cell in 
the accretion flow and the mean neutrino energy, respectively. For the neutrino emissivity 
distribution in the accretion flow, we adopt a simplified model: neutrinos and antineutrinos 
are emitted with a Fermi-Dirac distribution function with a zero-chemical potential and a 
temperature 



T cS (M,R) = T e s *(Mi 



lgn ; 




M < M ign 

M ign < M (57) 

(Mtrap,, c Mgn) 1/4 M>M tnw , e . 

where T*g(M,R) = (Q + /a) 1 ^. According to Zalamea & Beloborodov (2011), such a simpli- 
fication does not change the results significantly. Integrating l~Z v over the volume for 9 < 7r/4, 
where 9 is the angle from the rotation axis and R = — 30R g , we obtain the total energy 
deposition rate via neutrino pair annihilation above the accretion flow. 

We contrast the jet luminosity expected from the neutrino pair annihilation scenario 
with the jet luminosity generated by the BZ mechanism in Fig. [6] (thin lines). In all but 
the highest accretion rate state, the one in which all three neutrino flavors are very optically 
thick, we see that the BZ jet luminosity (for our fiducial parameters) is many orders of 
magnitude greater than the jet luminosity that would be produced by neutrinos. Even in 
that state, the neutrino luminosity is generally two orders of magnitude smaller than the BZ 
jet luminosity. In the accretion rate range most favorably to the BZ mechanism, the ratio is 
at least another order of magnitude larger. Thus, even if /(o/Mbh) is as small as ~ 10~ 2 , 
the BZ mechanism dominates at all accretion rates; in its favored range, it can be as small 
as ~ 10~ 3 and still generate a larger luminosity than the neutrino mechanism. 

When comparing the jet peak lu minosity with those observed in typi cal GRBs, ~ 
10 50 erg s _1 , after beaming corrections (IGhirlanda. Ghisellini fc Firmanill2006l ) we find that 



with our fiducial parameters those require accretion rate larger than 0.003 — O.OlM s 1 . The 
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Fig. 7. — Jet efficiency expected from the BZ mechanism. The parameters are the same as 
those in Fig. El 
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drop by a factor of ~ 5 in the efficiency at the transition Mi gn may be a mechanism that 
strongly modulates the jet power, leading to the observed variability in GRBs. 



4. Conclusion and Summary 

We have examined the Blandford-Znajek jet luminosity and efficiency expected from 
hyperaccreting black ho les, investigating the possibi lity that such a jet can account for ob- 



? 

served GRBs. Following iKrolik &: Piranl ( 120111 . 120121 ) . we estimated the Poynting luminosity 



using a simple dimensional argument: we assumed that the intensity of the magnetic field 
responsible for launching the jet is linked to the pressure in the innermost region of the 
accretion flow. By using an NDAF solution at the innermost region of the accretion flow, we 
have demonstrated that the luminosity of this Poynting-dominated jet can easily exceed the 
energy deposition rate expected from neutrino pair annihilation above the accretion flow. 
This means that MHD processes around the black hole horizon are a plausible mechanism 
for the formation of the relativistic jets required from the observations of GRBs. In ad- 
dition, we have derived the jet luminosity as a function of mass accretion rate and have 
shown that it has a step function-like behavior at a mass accretion rate corresponding to 
the transition between an advection-dominated accretion flow and a neutrino-dominated ac- 
cretion flow. This means that the jet power suddenly drops when the mass accretion rate 
decreases and crosses M ign . In the context of GRBs, this jump may correspond to switch- 
ing on and off the activity needed for internal shocks. It may also explain the steep decay 



a ^ 3 — 5; a = \d\ogF u / d\ogt\) often observe d after the prompt emission (ITagliaferri et al. 



20051 : iNousek et all 120061 : lO'Brien et all 120061 ) . We have also shown that the MHD jet for- 



mation process is most efficient when the accretion flow efficiently cools via optically-thin 
neutrino emission. In this regime, the jet luminosity is proportional to the mass accretion 
rate M, i.e., the efficiency does not depend on M. 

We stress that the BZ jet luminosity is larger than the energy deposition rate via neutrino 
pair annihilation whenever /(o/Mbh) <i 0.01. In the optimal accretion range of accretion 
rates, this lower bound is relaxed by at least another order of magnitude. Thus, for the BZ 
mechanism to dominate the neutrion mechanism, the central black hole must spin at least 
moderately fast, and the magnetic field geometry must be reasonably well-organized in the 
poloidal sense. However, as this lower bound on /(o/Mbh) shows, these thresholds are far 
from extreme. 
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